Dynamical phases and hysteresis in a simple one-lane traffic model 
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A two parameter model for single lane car-following is introduced and its equilibrium and non- 
equilibrium properties are studied. Despite its simplicity, this model exhibits a rich phenomenology, 
analogous to that observed in real traffic, like transitions between different dynamical regimes and 
hysteresis in the fundamental flux-density diagram. We show that traffic jams can spontaneously 
appear in clustered-like structures. In the jammed phase, we observe a slow relaxation phenomenon 
ruled by the outgoing car flux that determines the hysteretic dependence of the fundamental flux- 
density diagram. Coexisting phase regimes are also evidenced so as propagating or stationary density 
waves. The model can be easily calibrated to reproduce experimental observations. 

PACS numbers: 89.40. -f-k, 05.40. -f-j, 64.60. Cn, 64.60. Ht 
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The occurrence of traffic jams or small traffic conges- 
tions without obvious reasons are common effects that 
almost every driver has once experienced. These effects 
are typical signatures of the complex behavior of traffic 
flows and their study is of practical interest in the con- 
text of traffic control. The evolution of such spontaneous 
time-space structures has long attracted attention in the 
understanding of non-equilibrium properties of externally 
driven many-body systems. Sophisticated concepts were 
in particular developed to study critical phenomena and 
successfully applied to the description of phase transi- 
tions in traffic flow models 

Real traffic exhibits a very rich variety of phenomena. 
For this reason, the most adapted techniques for efficient 
traffic control are still debated. Experimental investiga- 
tions on highways revealed that traffic can exhibit well 
identified dynamical regimes that depend on the external 
car flux and on the car density |^||]. Phase transitions 
occur in traffic flows when the vehicle density exceeds a 
critical threshold. Below this threshold, traffic is free. 
Beyond it, vehicles either briefly slow down due to high 
density traffic or stop in a jam. Jams can appear with- 
out obvious reasons, they can merge and extend to large 
scales. The outcome of experimental measurements were 
moreover shown to be strongly influenced by the traf- 
fic behavior near the measurement site. Near ramps for 
example, peculiar behaviors such as avalanches or oscil- 
lations are known to take place ||]. Traffic flow results 
are thus far from universal 0| . 

Planning and optimizing real traffic flow has motivated 
many theoretical approaches, ranging from cellular au- 
tomata mode ls M to coupled maps M and from hydrody- 
namics §,|lO|-P| to kinetic theory Theoretical 
results depend on the drivers behavior. In deterministic 
models for instance non-linearities induce the jamming 
transition and the jammed phase persists in time p^ ] 
while for stochastic models a jamming transition occurs 



due to the intrinsic noise and no jam persists for ever 
Ip^l . Most of the theoretical models are based on heuris- 
tic arguments and their parameters have to be calibrated 
using experimental outputs. 

With the aim of retrieving experimental properties of 
traffic, we focus on a simple 2-parameters car following 
model. We consider a one lane street of length L — 10 
km. Vehicles on this street are all identical with size 
dear = 5 meters [T^ and indexed such that car num- 
ber i + 1 precedes car i. They cannot overtake and, due 
to periodic boundary conditions the vehicle ahead car 1 
is N. When approaching a slower car the driver slows 
down. The velocity of the front car remains unaffected. 
Car i has position Xi{t) and velocity Vi{t) < Vmax at time 
t, where Vmax is the maximal legal velocity. We assume 
that each driver determines his velocity in function of the 
headway to the front car. This is based on safety require- 
ments and on a simple empirical rule. Traffic security 
experts indeed often teach to beginners that driving con- 
ditions are safe at velocity Vi if the distance di with the 
vehicle ahead is such that di > Vi, where di is estimated 
in meters and Vi expressed in km/h. This means that at 
velocity 100 km/h a minimal headway security distance 
of 100 m is required. When following this rule, drivers 
have T = 3.6 seconds to react and avoid collisions. This 
time is independent of the velocity |l8|| and is supposed 
to include the front car stimulus, braking distance and 
drivers concentration level [p^ . The velocity law for car 
number i at time t is given by : 



Vmax for diit) > 25dcar (1) 
dt{t) for dear < di(t) < 25dcar 



where di{t) = mod{xi+i(t) — Xi{t),L) and Vi{t) is ex- 
pressed in km/h [^o| . d^ar and Vmax are the two param- 
eters of this model. For the latter to be consistent with 
experimental observations, we set Vmax = 125 km/h 
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The prefactor 25 in the first equation is then a conse- 
quence of the empirical driving rule. Indeed, in order 
to proceed safely, a driver at velocity Vmax has to slow 
down when the headway distance becomes smaller than 
125 m = 25d car- Although very simple, we show herein 
that simulations of Eqs. reproduce spatio-temporal 
observations of experimental traffic flows. For the time 
evolution of the position of each car we use a synchronous 
updating based on the Euler scheme Xi {t + At) = Xi (t) + 
Vi{t)At + 0{A€f with time step At = 0.1 s. Due to this 
time discretization and in order to avoid artificial car 
crash we have to constraint numerically v{i) = when 
d{i) < dear- This introduces a small dynamical noise 
that slightly perturbs our simulations. We checked how- 
ever that the main results presented herein still hold for 
At = 0.01 s. 

We study model (|^) under driven external car flux, ie. 
we allow the average car density on the street p to vary 
in time due to in-going (resp. out-going) cars through an 
on-ramp (resp. off-ramp) located at position x — km 
(resp. X = L = 10 km). To check the validity of this 
model, we first pay attention to the experimental obser- 
vations of reference For given position x on the street, 
we consider the one minute averaged velocity Vx of the 
cars crossing x and px the one minute averaged density in 
a 1 km band centered on x. We deal with one realization 
of the system and we model the external car flux by a sim- 
ple superposition of three Maxwellian functions. For this, 
we constraint the time evolution of the car density on the 
whole street to be p{t) = po J2i=i 2 3 6xp[— ((t — ti)/2)^], 
with ^1=7:00, t2=12:00 and t3 = 17:bo and where t is ex- 
pressed in hours. Using the experimental results of |^] 
we set pq ~ 50 cars/km and t^'s are supposed to be rush 
periods. Practically, we estimate the density p{t) every 
10 seconds and, according to its value, either keep the 
total number of cars to the same value or add/remove a 
single car from the street. This means that the maximal 
car flow is 360 cars/h which is consistent with constant 
external flux models |6| but below experimental observa- 
tions 1^. 

In Fig. ^ we display the local flux-density [q^, Px) fun- 
damental diagram and the time evolution of velocity ■ 
We consider two cross sections : The first is situated at 
position X = 5 km and the second near the off-ramp at 
X — 9 km. The {qx,Px) representation of Figs. || look 
very similar to the one obtained in Q. The numerical 
points tend however to accumulate in the neighborhood 
of (0 = 30 cars/km. This effect is a consequence of the 
simplicity of the time dependence of p{t) that tends to 
leave the system longer near p ~ 25 cars/km. It can be 
removed by adding small random fluctuations to p{t). 

When comparing W5 and vg, we observe larger fluctu- 
ations in the time behavior of vg (see full line in Fig. 
^(b) and (d)). They appear in time periods where the 
variations of p{t) are fast {i.e. when the external flux 
is large). These fluctuations are of particular interest 



near on-ramps where they can be seen as signature of 
avalanche-like effects that take place when cars arrive in 
a yet "busy" highway. 
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FIG. 1. Flux-density {qx,px) diagram (upper graphs) and 
average speed Vx (lower graphs full line) of the cars at crossing 
sections a; = 5 km (left column) and a; = 9 km (right column) 
during one minute intervals and for 24 hours. The dashed 
line in (b) and (d) corresponds to the time dependence of the 
total car density p{t) on the street. 

The effect of ramps was investigated both experimentally 
Q and theoretically [^,0. This local breakdown effect 
triggers jams and remains even if fluctuations in the ex- 
ternal flow are negligible. Time evolution of vi and qi is 
shown in Fig. |2| between time 4:00 and 13:00. Both fig- 
ures exhibit five different dynamical regimes: stable free 
traffic {SFT for iG [4:00,5:00]), local break down traf- 
fic {LBDT for te [5:00,6:30]), clustered traffic (CT for 
te[6:30,8:00]), synchronized traffic {ST for te[8:00,9:00]) 
and unstable free traffic {UFT for t€ [9:00,11:30]). In 
order to identify these regimes, we applied a noise to p{t) 
and observed that all the regimes but UFT were not 
significantely modified. UFT turns out to be sensitive 
to perturbations and disapears for too large noise am- 
plitudes jl3|. LBDT on the other hand appears when 
the density grows beyond p^i ~ 15 cars/km and is chara- 
terised by an avalanche- like deacreasing velocity vi for in- 
creasing density ]|ll| . This regime is located in the vicin- 
ity of the ramps and triggers complex non-homogeneous 
traffic structures that evolve over the whole street (see 
Fig. ^). CT consists in small jammed islands that appear 
away from ramps, persist in time and can be compared to 
droplets that form in supercooled gas. ST phase appears 
in highway measurements when a similar average car ve- 
locity is observed on all the lanes of the highway and is 
mostly localized in the neighborhood of the ramps. One 
line descriptions were however successfully used to study 
this regime J^0. ST is characterized by a low average 
velocity but a rather high value of the local car flux qi 
when compared to the one of the jammed traffic phase 
{JT). This latter phase does not appear in the previous 
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figures since po is too small. JT is a stable phase (any 
perturbation fade as time evolves) and is characterized 
by densely packed queues of cars that can extend in a 
single jam over the whole street. 
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FIG. 2. Detail of the time dependence of the average veloc- 
ity Hi (a), flux q\ and density pi (b) at cross section x = 1 km. 
Legends and simulation conditions are the same than in Fig. 
|l| The vertical dot-dashed line in (a) indicates the moment 
where the phase transition between free traffic flow and con- 
gested flow with critical density Pd ~ 15 cars/km. The dot- 
ted line in (a) corresponds to the time dependence of pi.In 
(b) we show the succession between SFT, LBDT, CT, ST 
and UFT. 

Many features of traffic flows can be compared to fluid 
dynamics. For instance JT can be seen as a compress- 
ible liquid without coexisting vapor whereas SFT can 
be trivially compared to a dilute gas Consequently, 
a criterion for model to be realistic constraints the 
transtion between JT and SFT to be histeretic as the 
liquid-gas transition. This prompts us to rapidly focus 
on the ensemble averaged total flux q and total density p. 
For each realization of the system, we perform a simula- 
tion that starts with p = 1 car/km, with a zero out flow 
4'out = of cars and a constant in flow 0i„ = ((> (in- flow 
period). Cars are introduced on the street one by one, 
to keep the system as close as possible to its stationary 
state, and q is averaged over the time running between 
two introductions. When density Pmax ^ 160 cars/km is 
reached, the street is uniformly jammed. At this density, 
in-flow is interrupted and we start to remove cars one by 
one setting 0i„ = and (j)out = 4'- We proceed this way 
down to density p = 1 car/km (out- flow period). This 
procedure is repeated over 100 realizations. Fig. || dis- 
plays the ensemble averaged (g, p) fundamental diagram 
for several values of cj). The SFT and the UFT regimes 
clearly show up for p < pc^ and p e [Pa, PC2] respectively, 
with pc2 ~ 20 cars/km. Both regimes do not strongly de- 
pend on the value of 4>. For p > Pd we also retrieve the 
characterisitic hysteresis loop of traffic behavior. The 
lower branch corresponds to in-flow periods. It is inde- 
pendent of the in-flow rate (j)i„. This is clearly not the 
case for out-flow periods as shown by the existence of 
several upper branches in Fig. 0. A direct consequence 



is that JT occurs beyond a threshold density that 
itself depends on the out-flow rate (j>out, where pc^ is de- 
fined as the density p > pc^ for which the two branches 
of the histeresis loop differ less than 10%. Hence, the 
global stability of model (|l]) is only determined by the 
out-going car flow (j>out ■ This property remains true even 
when modifying the model ]2l] ]. 
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FIG. 3. Flux-density fundamental diagram averaged over 
100 realizations for several values of the in/out car flow 0. 
Pc2 ~ 20 cars/km determines the threshold beyond which 
UFT appears. The lower branch (resp. upper branches) of 
the hysteresis loop is (resp. are) generated during the in-ffow 
(resp. out-flow) period, pc^ ~ 40, pc^ ~ 70 and pc^ ~ 110 
cars/km are deflned in the text and refer to simulations with 
(l> = 360 cars/h (full lines). For this value of (fi congested 
traffic appears when p £ [pc2,pcs] whereas JT appears for 
P > cars/km. 

For SFT, UFT and JT regimes the corresponding spa- 
tial distribution of the cars is uniform over the whole 
street. Since we focus on ensemble averaged quantities, 
only limited informations can be obtained from Figj^ on 
the complex space-time patterns of congested traffic. It is 
however interesting to note that for p > pc^ ~ 40 cars/km 
the in-flow branch is almost linear. This behavior is sim- 
ilar to the coexisting regime in fluid dynamics with the 
fluid and its gas in thermal equilibrium. Fig. |^ suggests 
that such a picture also holds for traffic flows although 
these systems are far from equilibrium. The explanation 
is that in this density range the system decomposes in 
two phases: A low density one [pcs < p < Pc^) that we 
identify as the CT and a high density one {p > pc^) cor- 
responding to JT. This is illustrated in Fig. ^ where 
we show the space-time density contour plot. Clustered 
structures of CT show up in Fig. § when pc^ < p < pc^ . 
As the density inside clusters is close to the critical value 
Pes this situation can be compared to the small droplets 
that form in a gas near condensation. As new cars arrive 
on the road, the JT first appears at time t=l:7 near the 
ramp and coexists with CT. As time evolves, due to in- 
creasing p, JT extends backwards and finally overcomes 
CT (at time 3:00 where p ~ Pc^)- Forward propagating 
density waves are evident in Fig. |4| The contour plot 
for p > pc^ in Fig.^ shows complex space time structure 
and an almost uniform JT phase that extends over the 
complete street when p ^ Pmax- 

Nucleated structures of CT also show up in the out- 
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flow period when is [6:7] {pc^ < p < pc^). But in this 
time range, it coexists with two quiet well separated JT 
and SFT phases. For longer times, when pc^ < p < Pa , a 
new separated phase between JT and SFT appears whit 
two density fronts. As p decreases with time, the in-flow 
of cars is no more sufficient to maintain JT that progres- 
sively disapears. In the down stream direction the front 
is standing while in the upstream one, the front propa- 
gates due to the 'evaporation' of the JT phase. From 
these results, Pa turns out to play the role of a physi- 
cal threshold between a triphasic JT-CT-SFT situation 
and a situation where JT and SFT coexist in a simple 
space-time shrinking structure. 
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FIG. 4. Time evolution of the total density p (left) and 
space-time density contour plot (right) for the run of Fig. |^ 
with (f) = 360 cars/h, during in flow t <4:30 and out-flow 
t >4:30 periods. In these flgures we restricted ourselves to 
the time period where p > pc^- In the left graph, the hor- 
izontal lines indicate the time at which the critical density 
Pa, i ~ 3, 4, 5, is reached. The street is divided in 20 sections 
of 500 meters and in each of them the density p is evaluated 
every minute. The legend of colors are given in cars/km. 

We considered herein a simple 2-parameters micro- 
scopic car-following model. When driven far from equi- 
librium by an external car flux, this model reproduces 
the main characteristics of experimental fundamental di- 
agram. Near ramps, we retrieve avanlanche-like effects 
so as synchronized traffic. The ensemble averaged funda- 
mental diagram exhibits the traditional histeresis loop of 
traffic models. This a signature of coexisting phases and 
we showed that among these phases two regimes have to 
be distinguished: a clustered phase where the system de- 
velops complex space time patterns and a regime where 
phases are separated by well identified density fronts. 
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